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We analysis the sea quark densities, based on the constituent quark model. To perform a direct fit with available experi- 
mental data, the parameterized Inverse Mellin technique is used. The calculation is extended to the NLO approximation 
for the singlet and non-singlet cases in DIS phenomena. We employ the approach of complete RG improvement(CORGI) 
where one is forced to identify and resum to all-orders RG-predictable ultraviolet logarithm terms which truly build the 
Q-dcpendence of QCD observable. The results are compared with the standard approach of perturbativc QCD in the 
MS scheme with a physical choice of scale. The results in the CORGI approach indicate a better agreement to the data. 



1. Introduction 

One of the feature of strong interaction field which 
is investigated by many phenomenological model is the 
Quark parton model. To construct the hadron structure 
from the parton densities, we use from the constituent 
quark model which seems to give us a better insight. 
A constituent parton is defined as a cluster of valence 
quarks accompanied by a cloud of sea quarks and 
gluons. They have been referred to as valons. It can 
be considered as a bound state in which, for instance, 
a proton consists of three valons, two U- valons and one 
D-valon which, on the one hand, interact with each other 
in a way that is characterized by the valon wave function 
and which on the other hand contribute independently 
in an inclusive hard collision with a Q 2 dependence that 
can be calculated in QCD at high Q 2 . These valons thus 
carry the quantum numbers of the respective valence 
quarks. Hwa ;2] found evidence for the valons in the 
deep inelastic neutrino scattering data, suggested their 
existence and applied it to a variety of phenomena. Hwa 
has also successfully formulated a treatment of the 
1ow-Pt reactions based on a structural analysis of the 
valons. Some papers can be found in which the valon 
model has been used to extract new information for 
parton distributions and hadron structure functions in 
unpolarized and polarized cases 

To improve and increase the reliability of the pertur- 
bative QCD calculation, we try to use the Complete 
Renormalization Group Improvement (CORGI) ap- 
proach which is related to the long standing problem of 
renormalization (/i)and factorization (M) scale depen- 
dence in QCD predictions. Usually in Renormalization 
Group(RG)-improving perturbation theory, it is assumed 
that these scales are related to the physical scale Q and 
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//= M=Q is normally chosen. The resulting fixed-order 
predictions depend on the choice of scales. If instead one 
insists that these dimentionful scales are independent 
of Q, one is forced to identify and resum to all-orders 
RG-predictable ultraviolet logarithms of Q which truly 
build the Q-dependcncc. In so doing all dependence on 
fi and M disappear. This Complete RG Improvement 
(CORGI) approach has previously been applied to the 
single scale case where there is only /z-dependence 0- 
It is then extended to the more complicated pattern of 
logarithms involved in the two scales problem of moments 
of structure functions [6j . In this approach the standard 
perturbative series of QCD observable is reconstructed 
in terms of scheme-invariant quantities. So it is expected 
to get more accurate results compared to those obtained 
in the standard perturbative QCD approach. 

Finally, in order to obtain directly all unknown parame- 
ters of the model, just by using the available experimental 
data, we use from the Inverse Mellin Technique, not in a 
numerical form but in a parameterized form. Usually one 
uses numerical computation to extract parton densities 
from their moments, but we believe that the parameter- 
ized Inverse Mellin technique is a mathematically power- 
ful tool which produces more reliable results than those 
obtained from the customary numerical method. Using 
the symmetry properties of the inverse Mellin transforma- 
tion and also a Taylor expansion of the integrand function, 
we are able to obtain the parameterized solution. 



2. Phenomenological constituent quark 
model 

The idea of quark cluster is not new. In this model, 
the hadron is envisaged as a bound state of valence 
quark clusters. For example the bound state of ir~ 
consists of a"anti-up" and "down" constituent quarks. 
In the static problems there is a little difference between 
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the usual constituent quarks and the valons, since the 
point-like nature of the constituent quarks is not a 
crucial aspect of the description, and has been assumed 
mainly for simplicity. But, in the scattering problems 
it is important to recognize that the valons, being 
clusters of partons, can not easily undergo scattering 
as a whole. The fact that the bound-state problem of 
the nucleon can be well described by three constituent 
quarks implies that the spatial extensions of the valons 
do not overlap appreciably. A physical picture of the 
nucleon in terms of three valons is then quite analogous to 
the usual picture of the deuteron in terms of two nucleons. 

To facilitated the phenomenological analysis the follow- 
ing simple form for the exclusive constituent quark inside 
the proton is assumed Q 

GuuD/ P {yi,y2,yz) = g(ym) a yi S(yi + 2/2 + 2/3 - 1), 

(i) 

where a and (3 are two free parameters and yi is the mo- 
mentum fraction of the i'th constituent quark. The U 
and D type inclusive constituent quark distributions can 
be obtained by double integration over the specified vari- 
ables: 



G u/ P (y) = J dy 2 J dy 3 G UUD/p (y,y 2 ,y 3 ) = 

gB(a+l,(3+l)y a (l-y) a +P+ 1 , (2) 



reasonable provided that the exchange of very soft gluons 
is responsible for the binding. Eq. (JjJJ involves also 
the assumption that in the deep inelastic scattering 
at high Q 2 the valons are independently probed, since 
the shortness of interaction time makes it reasonable 
to ignore the response of the spectator valons. Thus, 
through Eq. J^} we have broken up the hadron struc- 
ture problem into two parts. One part represented by 
G v / p {y), describes the wave functions of the proton in 
the valon representation. It contains all the hadronic 
complications due to the confinement. It is independent 
of Q 2 or the probe. The other part represented by 
F$(z = ^,Q 2 ), describes the virtual QCD processes of 
the gluon emissions and quark-pair creation. It refers 
to an individual valon independent of the other valons 
in the proton and consequently also independent of the 
confinement problem. It depends on Q 2 and the nature 
of the probe. 

Since the calculation in moment n-space is easier than 
the calculation in x-space, we work with to the moments 
of the distribution, defining 

M 2 , 3 (n,Q 2 ) = J o dxx n - 2 | ^ | (x,Q 2 ), (6) 

It then follows from Eq. (JSJ that 

M N (n,Q 2 ) = Y / M v/N (n)M v (n,Q 2 ). (7) 



G D / P (y) = / dyi J dy 2 G UUD / p {y 1 ,y 2 ,y) = 

gB(a+l,a + p + 2)yP(l-y) 2a+1 . (3) 

The normalization parameter g has been fixed by requir- 
ing 

Gu/p{y) dy= f G D/p {y) dy = 1, (4) 
o Jo 

and is equal to g = [B(a + 1, (3+ 1) B(a + 1, a + /3 + 2)]~\ 
where B(m, n) is the Euler-beta function. Hwa and Yang 
have recalculated the unpolarized valon distribution 
inside the proton with minimization method, using direct 
fit to CTEQ parton distributions, and have reported a 
new values of a and /3 which are found to be a = 1.76 
and p = 1.05. 

This model suggests that the structure function of a 
hadron involves a convolution of two distributions. Con- 
stituent quark distributions in the proton and structure 
function for each constituent quark so as 

F 2 p (x,Q 2 ) = Y / I dyG v/p (y)F^(z^- : Q 2 ). (5) 

„ Jx y 

We shall also assume that the three valons carry all 
the momentum of the proton. This assumption is 



The distributions we shall calculate (all referring to the 
proton) are those for sea quarks xu sea {x) = xu(x) — 
xd sea = xd(x) = xs(x) = xs(x) = • • •, which we shall 
generally denote by xq{x). Its moment is denoted by 
M sea (n,Q 2 ) or equivalently M sea (n, s) where s is evolu- 
tion parameter which in leading order is defined by 



s = log 



and generally at any required order is defined by 

<Ql) 



s = log 



a{Q 2 ) ' 



(8) 



(9) 



where a(Q 2 ) is related to the strong coupling constant by 
a(Q 2 ) = a s (Q 2 )/4:TT . Using Eq. 0, we obtain 

M sea (n,s) = —[2U(n)+D(n)] 
2/ 

x[M s (n,s)-M NS (n,s)], (10) 

where M° and M NS are singlet and non-singlet evolution 
function given in leading order by 

M NS (n,s) = exp(-d NS s), 



(11) 



1 



1 



M b (n,s) = -(l + p) cxp(-d+ s) + -(l-p) exp(-rf_.s 



(12) 
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These moments are the leading order solution of the 
renormalization group equation in QCD. The anomalous 
dimensions, d and other associated parameters are 
defined in |3j. 

We shall from now indicate the moment of sea quark 
densities by M(n,Q 2 ). In the NLO approximation, we 
haveH 



M Nb {n,Q 2 ) 



hi 

4tt ^ 2/3 2/3^ 

7 M /2/3o 



and 



M s (n,Q 2 ) 
M9(n,Q 2 ) 



P_ - 



1 a s (Q 2 ) - a s (Q 2 ) 



(13) 



P_. 7 .P_ 



2/3 4tt 
{QD *s(Q 2 ) (a s {Q 2 )\ {X+ - X - )/2Po 



-in 



4tt 



2fo + A+ - A_ 



a s (Q 2 ) 



A-/2/3 



(+ 



(14) 

The parameters (Vs, A's, /3's and P's) are defined in 
references 0, S E3, El E3| 

We will discuss how we can extend and improve the pre- 
cision of the calculation by direct use of available experi- 
mental data [13J , without referring to the input scales Qo 
and Ajjg. In this sense, we extract parton densities from 
the analytical moments of densities. A method devised 
to deal with this situation is to perform the integration of 
the inverse Mellin transformation in parameterized form 
which will be explained in section 4. 



3. Renormalization and factorization scale 
dependence 

The problem of renormalization scheme dependence in 
QCD perturbation theory remains on obstacle to making 
precise tests of the theory. It was pointed out [j| that the 
renormalization scale dependence of dimensionless physi- 
cal QCD observables, depending on a single energy scale 
Q, can be avoided provided that all ultraviolet logarithms 
which build the physical energy dependence on Q are 
resummed. This was termed complete Renormalization 
Group (RG)-improvement approach. For a single case of 
a dimensional observable 1Z(Q) with 



1Z(Q) =a + ria 2 + r 2 a 3 + ... + r n a n+1 + 



(15) 



The RS can be labelled by the non-universal coefficients 
of the beta- function and Ars fl^ . Self-consistency of 
perturbation theory that is the derivative of N-th order 
approximant 7Z(Q) with respect to the scheme labelled 
parameters, is of higher order than the approximant itself, 
will yields partial differential equations for coefficients r 2 , 
r3 and ... with respect to non-universal coefficients of 0- 
function and r\ coefficient. On integration of these partial 
differential equations, one finds 

V2{ri, c 2 ) = r 2 + cri + X 2 - c 2 



r 3 (ri,c 2 ,c 3 ) = n 3 + -cri 2 + {3X 2 - 2c 2 )ri + 



X 3 --c 3 



(16) 



General Structure is as follows 

r n {r 1 ,c 2 ,...,c„) = f n (r 1 ,c 2 , . . .,c n _i) + X n 
-c n /(n - 1); 

where X n are Q-independent and RS-invariant and are 
unknown unless a complete N n LO calculation has been 
performed. Now we can reformulated R{Q) as it follows 

K(Q) = a + r ia 2 + (r 2 + cr x + X 2 - c 2 )a 3 + ... . (17) 

Given a NLO calculation, r\ is known but X 2 , X3, . . . are 
unknown. Thus the complete subset of known terms in 
Eq. (JUJ at NLO is 

a =a + no 2 + (rj + en - c 2 )a 3 + ... , (18) 

and it is RS-invariant. Choose r\ — 0, c 2 — C3 = = 
c n = we obtain a = a 6]. At NNLO calculation X 2 is 
unknown. Further infinite subset of terms are known and 
can be resummed to all orders, 



AW = X 2 a 3 + 3A 2 na 4 + . . . . 
Finally we will arrive at 
11{Q) = a + A 2 a 3 + A 3 ao 4 + - • . + A„a " +1 - 



(19) 



(20) 



which in fact is the expansion of QCD observable in 
Complete RG-Improvement (CORGI) approach. Here 
ao = a(0, 0,0, ) is the coupling in this scheme and sat- 
isfies 



1 



, cln 

ao V 1 + ca o 



ca 



bin 



Q_ 

An 



(21) 



In fact the solution of this transcendental equation can 
be written in closed form in terms of the Lambert W- 
function 

[HD3, defined implicitly by W{z) exp(VF(z)) = 



1 



ao = 



z(Q) 



c[l + W(z(Q))] 
1 ( Q_ n - b/c 
An, 



(22) 
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where 6 and c are the first two universal terms of QCD 
/3-function. 

Each term in Eq. (|2U|) involves a resummation of 
infinite terms at specified order. For instance, the first 
term in this equation, ao, is representing a resummation 
over NLO contributions and the second term, A^ajj, a 
resummation over the NNLO contributions of all terms 
in Ea. l|15|l . The advantage of CORGI approach is that 
not only each term in the related perturbative series is 
scheme independent but also because it involves a re- 
summation of RG- predictable terms, it should yield more 
accurate approximations than the standard truncation of 
Eq. HS1). 

If we intend to employ the CORGI approach to 
extract sea quark densities in the LO approximation 
we need just to change the RG-coupling constant a 
to ao which exists in definition of evolution param- 
eter s that is defined by s — log °|gSj| ■ We should 
note that in the definition of ao, the A-r. quantity will 
depend on n and this coefficient is observable-dependent. 

To do the NLO calculation in the standard approach 
to extract sea quark densities, it is necessary to have 
anomalous dimensions and Wilson coefficient functions 
for the singlet and non-singlet sectors involved, for 
instance, deep inelastic of lepton-nucleon scattering 
or e + e~ annihilation. According to Eqs. I|13I14|) and 
following 0, U 0, [n| , the results for singlet and non- 
singlet moments can be obtained and the final results 
for sea quark density will be independent of the chosen 
observable. 

To employ the CORGI approach in higher order, we 
should use the general form for moments |l7| : 

M(n,Q 2 ) = 

/ coo^ (l + X 2 (n)aUn) (23) 

VI +a {n)J 

+X 3 (n)a 3 (n) + ■■■ + X k {n)a k Q (n) + •••), 

where ao(n) is defined by Eq. (1221) and Xk(n) are the 
scheme invariant constants which where introduced be- 
fore. For the case of dependence on [i and M scales we 
have [13 

{2cY c/b fdAn) n(n)\ 

where di(n) is the MS NLO anomalous dimension co- 
efficient, and r\{n) is computed in the MS scheme with 
/i = Q. The (2c/6)~ C//fc factor corresponds to the standard 
convention for defining Xjjg. In the NLO approximation 
of the CORGI approach , we just need to keep the first 



two terms in the above series. The A^n) is defined by 



X 2 {n) 



■1 



2 2d(n 
c 2 d(n) 



■r\{n) 



cd\{n) 0*2 



26 



2b 



2b 



ri(n). 



(25) 



The NLO CORGI invariant X 2 (n) can be computed 
from the MS results for ri(n), ra (n) , aj (n) & d-2(n) 
[l0| . We should note that in Refs.pjL Il7| the adopted 
conventions for defining the anomalous dimensions and 
coefficient functions are different with respect to other 
references. To obtain M(n, Q 2 ) for the singlet case, we 
need to diagonalize the anomalous dimension matrix. To 
avoid this we restrict ourselves just to considering the 
non-singlet case. 



4. Fitting methods- Parameterized Solution 

In order to check the validity of the valon model and 
also to increase its ability to get more reliable parton 
densities (sea quark densities in our case), we make use 
of a method which we call /3-fitting. 

We assume the following phenomenological form for the 
sea quark densities: 



xq(x, Q 2 ) = a x b (1 - x) c (1 + dx + eyfx). 



(26) 



The parameters a, b, c, d and e are generally Q 2 - 
dependent. The motivation for choosing this functional 
form is that the term x b controls the low- a; behavior of 
sea densities, and (1 — x) c that at large values of x. The 
remaining polynomial factor accounts for the additional 
medium- :r values. 

The results of calculation indicates that the chosen 
functional form will yield a better fitting for the moment 
of the distribution, than that assumed for the sea quark 
distribution in Ref@|- 



Using the Mellin transformation 



M{n,Q 2 )= / x n ~ 2 xq(x,Q 2 )dx, 
Jo 



(27) 



where xq(x,Q 2 ) has been defined in Eq. I|26|l . we arrive 
at 



M(n,Q 2 ) = «r(l + c) 



r(-l + 6 + n) 
T(6 + c + n) 



er(-i + 6 + n) dT(b + n) 



r(i+6 + c + n) r(l + 6 + c + n) 



(28) 
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The unknown parameters a, b, c, d and e in Eq. I|28l) are 
obtained from the fitting of analytical results of moments 
in Eq. I|1(J|I over the parameterized form of this equation 
which is in terms of T or eventually Euler Beta functions. 
So we call this method /3-fitting. The quantities M s 
and M NS in Eq. 0, are known analytically from QCD 
calculation and moments of valon distributions, U(n) 
and D(n), are defined in 0. 

Now to report directly all exist parameters, just by us- 
ing the available experimental data, we can use from the 
Inverse Mellin Technique, not in a numerical form but in a 
parameterized form. For this propose, we take the Inverse 
Mellin of moment of distribution in complex space: 

1 fC+lOO 1 

F(x,Q 2 ) = — —M(n,Q*). (29) 

zm Jc—ioo x 

Since the variable c should be placed on the right hand 
side of all singularities and considering this point that 
all singularities of moments will occur for n less than 
2, so we choose c equal to 3. The integrated interval 
was first [3 — lOi, 3 + 10z]. If we choose the interval 
[3 — 20i,3 + 20i], there will be a difference only of order 
about 10~ 5 with respect to last interval. 

Since the function M(n, Q 2 ) is not a simple function, 
our machine will not be able to compute Eq. I|29|) . We 
use from this point that this integral with respect to real 
axis, is symmetric. So first we do integral for interval 
[c = c+Oi, c+mi] where the final result is twice this result. 

The technique which we used to integrate the M(n, Q 2 ) 
in 129|) is that we first choose a small interval, say [c + 
Oi, c + 2ei]. Then we expand M(n, Q 2 ) about c + ei. If 
e is small enough, we can keep just the first term of this 
expansion. In next step, we repeat the calculation for 
interval [c+2ei, c±4ei] and expand M(n, Q 2 ) about c+3ei. 
Repeating this procedure and adding all results, we are 
able to calculate Eq. (|29|l completely in parameterized 
form. 

Now we are in a position (by direct fitting of Fix, Q 2 ) 
from Eq. (|29|l over available experimental data [13) to 
obtain the unknown parameters of this function which 
include a and (3 (valon parameters) and Qo and A-^-j. 
These last two parameters occur in the definition of 
the evolution parameter, s. Unfortunately, from this 
direct fitting, we were not able to get reliable results for 
the unknown parameters. To overcome this difficulty, 
we inserted the values of the valon parameters a and 
f3, quoted from 0, in the distribution function, so 
this function will just depend on the Qo and A^j^ 
parameters, / = f(x,Qo,Ajjg). To get the best fitting 
X 2 value, we multiplied this function by an auxiliary 
term, Ax B (1 — x) c and fitted to obtain the unknown 
parameters Qo, A-p-5, A, B and C. We did the fitting 
for , Q 2 = 3 GeV 2 and Q 2 = 5 GeV 2 where we assumed 
a number of active quark flavours Nf — 3 and took an 



Table I Numerical values of the unknown parameters, 
resulted from the parameterized solution of the inverse 
Mellein transform technique in the LO approximation. 



Fitting parameters 


A 


3.358 ±0.124 


B 


0.134 ±0.033 


C 


-3.345 ± 0.264 


Qo 


0.422 ± 0.009(Gev) 




0.252 ± 0.004(Gev) 



average of the fitted parameters. The difference between 
the fitted parameters can be reported as an error of the 
calculations. The results are tabulated in Table I. 

One way to check the validity of calculations is to ex- 
tract the quark densities using the phenomcnological form 
in Ea. (|26(l . For this purpose, we went back to the /3-fitting 
method and repeated those calculations, but this time 
with the extracted values of Qo and A-^-<j (a and (3 taken 
from Ref 7]. As a consequence, new values for parameters 
contained in Ea. l2tj|) will be obtained. In order to get the 
best fit, we again multiply Ea. l2tj|) by the auxiliary term 
A x B (1 — x) c , as we did in the inverse Mellin Technique. 
We refer to this result as the "improved Inverse Mellin 
(IM) technique" . These new results in two different stan- 
dard and CORGI approaches and their comparison with 
experimental data are plotted in Fig.l and Fig. 2. As can 
be seen from Fig. 2 the result in CORGI approach will 
show better behavior in small x values. 



5. Conclusions 

Constituent quark model or valon model as a good 
candidate to describe deep inelastic scattering and to ex- 
tract sea quark densities inside the nucleon is used. The 
model bridges the gap between the bound state problem 
and the scattering problem for hadrons. It is possible 
to use this model in the LO and NLO approximations, 
using two different standard and CORGI approaches. 
The most important motivation for the CORGI approach 
is that by completely resumming all the UV logarithms, 
one correctly generates the physical dependence of the 
moments Ai(n, Q 2 ) on the DIS energy scale Q. 

In order to get the unknown parameters of the model 
from the fitting over the available experimental data, 
parameterized Inverse Mellin technique is used. The 
compared results will show that the CORGI approach 
indicates better consistency with available experimental 
data. An alternative "/3-fitting" method which is based 
on using the defining parameters of the Valon model, was 
also used to confirm the calculations of the parameterized 
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Q* = 3 GeV 2 
(LO appro.) 



Data from Ref.[13] 
CTEQ6 

CORGI-IM improved 
Standard-IM improved 




Q 2 = 5 GeV 2 
(LO appro.) 



Data from Ref.[13] 
CTEQ6 

CORGI-IM improved 
Standard-IM improved 




Figure 1: The results of Inverse Mellin (IM) technique for the standard and CORGI approaches in the LO approximation. 





Q 2 = 3 GeV 2 




(NLO appro.) 




□ Data from Ref.[13] 




CTEQ6 




IM-CORGI 


- 1 \t 


IM-Standard 











0.1 0.2 0.3 0.4 0.5 0.6 

x-Bjorken 





Q 2 = 5 GeV 2 




(NLO appro.) 




□ Data from Ref.p 3] 




CTEQ6 




IM-CORGI 




IM-Standard 










1 0.2 0.3 0.4 0.5 




x-Bjorken 



Figure 2: The results of Inverse Mellin (IM) technique for NLO approximation in both of Standard and CORGI approach 
Inverse Mellin technique. 



These calculations can be extended to the higher 
order of standard and CORGI approaches, using the 
recent analytical calculations for Wilson coefficients and 
anomalous dimensions. 
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We analysis the sea quark densities, based on the constituent quark model. To perform a direct fit with available 
experimental data, the parameterized Inverse Mellin transform technique is used. The calculation is extended to the 
NLO approximation for the singlet and non-singlet cases in DIS phenomena. We employ the approach of complete RG 
improvement(CORGI) where one is forced to identify and resum to all-orders RG-prcdictable ultraviolet logarithms of 
Q which truly build the Q-dependence. The results are compared with the standard approach of perturbative QCD in 
the MS scheme with a physical choice of scale. The results in the CORGI approach indicate a better agreement to the 
data. 



1. Introduction 

One of the feature of strong interaction field which 
is investigated by many phenomenological model is the 
Quark parton model. To construct the hadron structure 
from the parton densities, we use from the constituent 
quark model Q which seems to give us a better insight. 
A constituent parton is defined as a cluster of valence 
quarks accompanied by a cloud of sea quarks and 
gluons. They have been referred to as valons. It can 
be considered as a bound state in which, for instance, 
a proton consists of three valons, two U-valons and one 
D-valon which, on the one hand, interact with each other 
in a way that is characterized by the valon wave function 
and which on the other hand contribute independently 
in an inclusive hard collision with a Q 2 dependence that 
can be calculated in QCD at high Q 2 . These valons thus 
carry the quantum numbers of the respective valence 
quarks. Hwa found evidence for the valons in the 
deep inelastic neutrino scattering data, suggested their 
existence and applied it to a variety of phenomena. Hwa 
[2j has also successfully formulated a treatment of the 
low-Pr reactions based on a structural analysis of the 
valons. Some papers can be found in which the valon 
model has been used to extract new information for 
parton distributions and hadron structure functions in 
unpolarized and polarized cases 3 . 

To improve and increase the reliability of the pertur- 
bative QCD calculation, we try to use the Complete 
Renormalization Group Improvement (CORGI) ap- 
proach which is related to the long standing problem of 
renormalization (/i)and factorization (M) scale depen- 
dence in QCD predictions. Usually in Renormalization 
Group(RG)-improving perturbation theory, it is assumed 
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that these scales are related to the physical scale Q and 
/i= M=Q is normally chosen. The resulting fixed-order 
predictions depend on the choice of scales. If instead one 
insists that these dimcntionful scales are independent 
of Q, one is forced to identify and resum to all-orders 
RG-predictable ultraviolet logarithms of Q which truly 
build the Q-dependence. In so doing all dependence on 
(i and M disappear. This Complete RG Improvement 
(CORGI) approach has previously been applied to the 
single scale case where there is only /i-dependence 
It is then extended to the more complicated pattern of 
logarithms involved in the two scales problem of moments 
of structure functions . In this approach the standard 
perturbative series of QCD observable is reconstructed 
in terms of scheme-invariant quantities. So it is expected 
to get more accurate results compared to those obtained 
in the standard perturbative QCD approach. 

Finally, in order to obtain directly all unknown parame- 
ters of the model, just by using the available experimental 
data, we use from the Inverse Mellin Technique, not in a 
numerical form but in a parameterized form. Usually one 
uses numerical computation to extract parton densities 
from their moments, but we believe that the parameter- 
ized Inverse Mellin technique is a mathematically power- 
ful tool which produces more reliable results than those 
obtained from the customary numerical method. Using 
the symmetry properties of the inverse Mellin transforma- 
tion and also a Taylor expansion of the integrand function, 
we are able to obtain the parameterized solution. 



2. Phenomenological constituent quark 
model 

The idea of quark cluster is not new. In this model, 
the hadron is envisaged as a bound state of valence 
quark clusters. For example the bound state of ir~ 
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consists of a"anti-up" and "down" constituent quarks. 
In the static problems there is a little difference between 
the usual constituent quarks and the valons, since the 
point-like nature of the constituent quarks is not a 
crucial aspect of the description, and has been assumed 
mainly for simplicity. But, in the scattering problems 
it is important to recognize that the valons, being 
clusters of partons, can not easily undergo scattering 
as a whole. The fact that the bound-state problem of 
the nucleon can be well described by three constituent 
quarks implies that the spatial extensions of the valons 
do not overlap appreciably. A physical picture of the 
nucleon in terms of three valons is then quite analogous to 
the usual picture of the deuteron in terms of two nucleons. 

To facilitated the phenomenological analysis the follow- 
ing simple form for the exclusive constituent quark inside 
the proton is assumed Q 

GuuD/piVi, 2/2,2/3) = 9 (2/12/2)" 2/3 6 (yi +J/2 + 2/3 - 1), 

where a and (3 are two free parameters and yi is the mo- 
mentum fraction of the i'th constituent quark. The U 
and D type inclusive constituent quark distributions can 
be obtained by double integration over the specified vari- 
ables: 



Gu/ P {y) = J dy 2 J dy 3 G UUD/p (y,y 2 ,y 3 ) = 

gB(a+l,j3+l)y a (l-y) a+ P + \ (2) 



We shall also assume that the three valons carry all 
the momentum of the proton. This assumption is 
reasonable provided that the exchange of very soft gluons 
is responsible for the binding. Eq. (J5J involves also 
the assumption that in the deep inelastic scattering 
at high Q 2 the valons are independently probed, since 
the shortness of interaction time makes it reasonable 
to ignore the response of the spectator valons. Thus, 
through Eq. JSJ) we have broken up the hadron struc- 
ture problem into two parts. One part represented by 
G v i p (y), describes the wave functions of the proton in 
the valon representation. It contains all the hadronic 
complications due to the confinement. It is independent 
of Q 2 or the probe. The other part represented by 
■^2(2 = yiQ 2 )i describes the virtual QCD processes of 
the gluon emissions and quark-pair creation. It refers 
to an individual valon independent of the other valons 
in the proton and consequently also independent of the 
confinement problem. It depends on Q 2 and the nature 
of the probe. 

Since the calculation in moment n-space is easier than 
the calculation in x-space, we work with to the moments 
of the distribution, defining 

M 2 , 3 (n,Q 2 ) = J dxx n - 2 I ^ J (x,Q 2 ), (6) 

It then follows from Eq. JSJ that 

M N {n,Q 2 ) = Y. M v/N(n)M v (n,Q 2 ). (7) 



G D / P (y) = J dyi J dy 2 G UUD / p (y 1 ,y 2 ,y) = 

gB(a+l,a + p + 2)yf 3 (l-y) 2a+1 . (3) 
The normalization parameter g has been fixed by requir- 



ing 



Gu/ P {y) dy = G D/p (y) dy = 1, 



(4) 



and is equal to g = [B(a + 1 , f3+ 1) B(a + 1, a + /3 + 2)] _1 , 
where B(m, n) is the Euler-beta function. Hwa and Yang 
7] have recalculated the unpolarized valon distribution 
inside the proton with minimization method, using direct 
fit to CTEQ parton distributions, and have reported a 
new values of a and (3 which are found to be a = 1.76 
and (3 = 1.05. 

This model suggests that the structure function of a 
hadron involves a convolution of two distributions. Con- 
stituent quark distributions in the proton and structure 
function for each constituent quark so as 



Q 2 ) = E / d yGv/ P (y) F£(z = -,Q 2 
,, Jx y 



(5) 



The distributions we shall calculate (all referring to the 
proton) are those for sea quarks xu sea (x) = xu(x) — 
xd sea = xd(x) = xs(x) = xs(x) = ■ ■ •, which we shall 
generally denote by xq{x). Its moment is denoted by 
M sea (n,Q 2 ) or equivalently M sea (n, s) where s is evolu- 
tion parameter which in leading order is defined by 



and generally at any required order is defined by 

<Ql) 



S = loe 



a(Q 2 ) ' 



(8) 



(9) 



where a(Q 2 ) is related to the strong coupling constant by 
a(Q 2 ) = a s {Q 2 )/An . Using Eq. 0, we obtain 

M sea (n,s) = jj[2U(n) + D(n)} 

x [M s (n, s) - M NS {n, s)] , (10) 

where M° and M NS are singlet and non-singlet evolution 
function given in leading order by 



M NS {n,s) = exp{-d NS s), 



(11) 
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M s {n,s) = \{l+p) exp(-d+s) + i(l-p) exp(-d_«). 

(12) 



These moments are the leading order solution of the 
rcnormalization group equation in QCD. The anomalous 
dimensions, d and other associated parameters are 
defined in Q. 

We shall from now indicate the moment of sea quark 
densities by M(n,Q 2 ). In the NLO approximation, we 
haveS 



M NS (n,Q 2 ) = [ 1+ MQVMQ^) (ff - ffi) 



and 



(13) 



Group (RG)-improvement approach. For a single case of 
a dimensional observable 1Z(Q) with 



K(Q) 



r 2 a 



+ r n a 



n+l 



(15) 



The RS can be labelled by the non-universal coefficients 
of the beta- function and Ars [Til- Self-consistency of 
perturbation theory that is the derivative of N-th order 
approximant 1Z(Q) with respect to the scheme labelled 
parameters, is of higher order than the approximant itself, 
will yields partial differential equations for coefficients r 2l 
r3 and ... with respect to non-universal coefficients of 0- 
function and r\ coefficient. On integration of these partial 
differential equations, one hnds Q 



f*2(n,c 2 ) = r 2 + cri + X 2 



C2 



r3{n 1 c 2 ,c 3 ) = ri 3 + -cri 2 + (3A 2 ~ 2c 2 )n 
X,-\c, 



'M s (n,Q 2 )\_ 
MS( n ,Q 2 ) 



p 1 a s {Ql)-a s (Q 2 ) p _ r;p 



2/3o 



4tt 



MQo) a s {Q 2 ) ( a s (Q 2 ) 
^ WQo) 



-in 



P_. 7 .P 4 



(A + -A_)/2/V 



a s (Q 2 ) 



:2 ,, A_/2ft> 



•(+ 



")■ 

(14) 

The parameters (Vs, A's, /3's and P's) are defined in 
references @, 13 Efl El E3 ■ 

We will discuss how we can extend and improve the pre- 
cision of the calculation by direct use of available experi- 
mental data ^3 j without referring to the input scales Qo 
and Ajjg. In this sense, we extract parton densities from 
the analytical moments of densities. A method devised 
to deal with this situation is to perform the integration of 
the inverse Mellin transformation in parameterized form 
which will be explained in section 4. 



3. Renormalization and factorization scale 
dependence 

The problem of renormalization scheme dependence in 
QCD perturbation theory remains on obstacle to making 
precise tests of the theory. It was pointed out p| that the 
renormalization scale dependence of dimensionless physi- 
cal QCD observables, depending on a single energy scale 
Q, can be avoided provided that all ultraviolet logarithms 
which build the physical energy dependence on Q are 
resummed. This was termed complete Renormalization 



i . (16) 

General Structure is as follows 

r n (ri, c 2 , . . ., On) = f n (ri,C2,...,Cn-l)+X n 
-c n /(n - 1); 

where X n are Q-independent and RS-invariant and are 
unknown unless a complete N n LO calculation has been 
performed. Now we can reformulated R(Q) as it follows 



Tl(Q) =a + ria 2 + (r\ + en + X 2 - c 2 )a 3 + 



(17) 



Given a NLO calculation, r\ is known but X 2 , X 3l . . . are 
unknown. Thus the complete subset of known terms in 
Eq. (EJ at NLO is 



a n — a + r\a 2 + {r\ + cr\ — c 2 )a 3 + 



(18) 



and it is RS-invariant. Choose r\ — 0, c 2 — C3 = = 
c n = we obtain a — a 6]. At NNLO calculation X 2 is 
unknown. Further infinite subset of terms are known and 
can be resummed to all orders, 

A 2 a 3 = X 2 a 3 + 3A 2 na 4 + . . . . (19) 

Finally we will arrive at 

K{Q)=a Q +X 2 ao 3 + X 3 a Q 4 + ... + X n a n+1 + ... , (20) 

which in fact is the expansion of QCD observable in 
Complete RG-Improvement (CORGI) approach. Here 
a = <z(0,0, 0, ) is the coupling in this scheme and sat- 
isfies 



1 + C ln(-^- 

ao \ 1 + ca 



bin 



A 



■R 



(21) 
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In fact the solution of this transcendental equation can 
be written in closed form in terms of the Lambert dis- 
function defined implicitly by W(z) exp(W(z)) = 
z, 

1 

fl ° ~ ~ c[l + W(z(Q))} 
Z[ Q) = -\(£f hlC . (22) 

where b and c are the first two universal terms of QCD 
/3- function. 

Each term in Eq. (1201) involves a resummation of 
infinite terms at specified order. For instance, the first 
term in this equation, ao, is representing a resummation 
over NLO contributions and the second term, A^af], a 
resummation over the NNLO contributions of all terms 
in Ea. l|15|l . The advantage of CORGI approach is that 
not only each term in the related perturbative series is 
scheme independent but also because it involves a re- 
summation of RG-predictable terms, it should yield more 
accurate approximations than the standard truncation of 
Eq. GSJ. 

If we intend to employ the CORGI approach to 
extract sea quark densities in the LO approximation 
we need just to change the RG-coupling constant a 
to ao which exists in definition of evolution param- 
eter s that is defined by s = log ^^§y ■ We should 
note that in the definition of ao, the A-r. quantity will 
depend on n and this coefficient is observable-dependent. 

To do the NLO calculation in the standard approach 
to extract sea quark densities, it is necessary to have 
anomalous dimensions and Wilson coefficient functions 
for the singlet and non-singlet sectors involved, for 
instance, deep inelastic of lepton-nucleon scattering 
or e + e~ annihilation. According to Eqs. I|13I14|) and 
following [j| 0, 0, [nj , the results for singlet and non- 
singlet moments can be obtained and the final results 
for sea quark density will be independent of the chosen 
observable. 

To employ the CORGI approach in higher order, we 
should use the general form for moments |17| : 

M{n,Q 2 ) = 

(f \ \ d(n)/b 
t^P^t (l + X 2 (n)og(n) (23) 

l + a (n)J 

+X 3 (n)a 3 (n) + ■■■ + X k {n)a k Q {n) + ••■), 

where ao(n) is defined by Eq. 1|22[) and X^{n) are the 
scheme invariant constants which where introduced be- 
fore. For the case of dependence on /i and M scales we 



have [13 

/2cV c/h fdxin) n(n)\ , , 

where di(n) is the MS NLO anomalous dimension co- 
efficient, and ri(n) is computed in the MS scheme with 
/i = Q. The (2c/6)~ c / b factor corresponds to the standard 
convention for defining Xjj-g. In the NLO approximation 
of the CORGI approach , we just need to keep the first 
two terms in the above series. The X2(n) is defined by 

m 

/ . —1 b 9/ . cdi(n) 0*2 (n) 

X ^ n) = --2d^ ) rl{n) ~-k 1 + ^b L + 

C -^+r 2 (n). (25) 

The NLO CORGI invariant X^ri) can be computed 
from the MS results for ri(n), r 2 (n), di(n) & g?2(™) 
[Toj . We should note that in Refs.jE 113 tne adopted 
conventions for defining the anomalous dimensions and 
coefficient functions are different with respect to other 
references. To obtain M(n, Q 2 ) for the singlet case, we 
need to diagonalize the anomalous dimension matrix. To 
avoid this we restrict ourselves just to considering the 
non-singlet case. 



4. Fitting methods- Parameterized Solution 

In order to check the validity of the valon model and 
also to increase its ability to get more reliable parton 
densities (sea quark densities in our case), we make use 
of a method which we call /3-fitting. 

We assume the following phenomenological form for the 
sea quark densities: 

xq(x, Q 2 ) = ax b {l-x) c (l + dx + e^fx). (26) 

The parameters a, b, c, d and e are generally Q 2 - 
dependent. The motivation for choosing this functional 
form is that the term x b controls the low- a; behavior of 
sea densities, and (1 — x) c that at large values of x. The 
remaining polynomial factor accounts for the additional 
medium- a; values. 

The results of calculation indicates that the chosen 
functional form will yield a better fitting for the moment 
of the distribution, than that assumed for the sea quark 
distribution in Ref|2|. 

Using the Mellin transformation 

M(n,Q 2 )= [ x n ~ 2 xq(x,Q 2 )dx, (27) 
Jo 
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where xq(x, Q 2 ) has been defined in Eq. 
at 



est, 



we arrive 



M(n, Q 2 ) 



aF(l + c) 

, er(-i 



T(-l + b + n) 
T(b + c + n) 
b + n) dT(b + n) 



r(i + 6 + c + n) T(l + b + c + n) 



(28) 

The unknown parameters a, b, c, d and e in Eq. I|28|) are 
obtained from the fitting of analytical results of moments 
in Eq. 1)11) [I over the parameterized form of this equation 
which is in terms of T or eventually Euler Beta functions. 
So we call this method /3-fitting. The quantities M s 
and M NS in Eq. {T® , are known analytically from QCD 
calculation and moments of valon distributions, U(n) 
and Din), are defined in pj. 

Now to report directly all exist parameters, just by us- 
ing the available experimental data, we can use from the 
Inverse Mellin Technique, not in a numerical form but in a 
parameterized form. For this propose, we take the Inverse 
Mellin of moment of distribution in complex space: 



F(x,Q 2 ) = ^- I 
2m J c . 



c-\-ioo 



dn 



Y M(n, Q ). (29) 



Since the variable c should be placed on the right hand 
side of all singularities and considering this point that 
all singularities of moments will occur for n less than 
2, so we choose c equal to 3. The integrated interval 
was first [3 — KM, 3 + 10z]. If we choose the interval 
[3 — 20i,3 + 20i], there will be a difference only of order 
about 10~ 5 with respect to last interval. 

Since the function M(n,Q 2 ) is not a simple function, 
our machine will not be able to compute Eq. i|29|) . We 
use from this point that this integral with respect to real 
axis, is symmetric. So first we do integral for interval 
[c = c+Oi, c+mi] where the final result is twice this result. 

The technique which we used to integrate the 
M(n, Q 2 ) in (JUJl is that we first choose a small interval, 
say [c + Oi,c + 2ei}. Then we expand M(n,Q 2 ) about 
c + ei. If e is small enough, we can keep just the first 
term of this expansion. In next step, we repeat the 
calculation for interval [c + 2ei, c + 4ei] and expand 
M{n, Q 2 ) about c + 3ei. Repeating this procedure and 
adding all results, we are able to calculate Eq. i(2*9")l 
completely in parameterized form. 

Now we are in a position (by direct fitting of F(x. Q 2 ) 
from Eq. (|29|l over available experimental data [jjj) to 
obtain the unknown parameters of this function which 
include a and f3 (valon parameters) and Qq and A-p-^. 
These last two parameters occur in the definition of 
the evolution parameter, s. Unfortunately, from this 






Q* = 5 GeV 2 


-\ 


(LO appro.) 


■ 1 


DatetromRef.113] 


- V 

■ \ 




\ 
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Figure 1: The results of Inverse Mellin (IM) technique for the 
standard and CORGI approaches in the LO approximation. 




Q* = 5 GeV 1 
(NLO appro.) 



Data trom Ref.[1 3| 
CTEQ6 
IM-CORGI 
IM-Standard 




Figure 2: The results of Inverse Mellin (IM) technique for 
NLO approximation in both of Standard and CORGI 
approach. 



direct fitting, we were not able to get reliable results for 
the unknown parameters. To overcome this difficulty, 
we inserted the values of the valon parameters a and 
j3, quoted from 0, in the distribution function, so 
this function will just depend on the Qq and ^.jjg 
parameters, / = f(x,Qo,Aj^g). To get the best fitting 
X 2 value, we multiplied this function by an auxiliary 
term, Ax B (1 — x) c and fitted to obtain the unknown 
parameters Qo, ^ms^ ^ ^ an< ^ ^ ■ ^ e ^ e fitting 
for , Q 2 — 3 GeV 2 and Q 2 = 5 GeV 2 where we assumed 
a number of active quark flavours Nf — 3 and took an 
average of the fitted parameters. The difference between 
the fitted parameters can be reported as an error of the 
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Table I Numerical values of the unknown parameters, 
resulted from the parameterized solution of the inverse 
Mellein transform technique in the LO approximation. 



Fitting parameters 



A 3.358 ±0.124 

B 0.134 ±0.033 

C -3.345 ± 0.264 

Qo 0.422 ± 0.009(Gev) 

A 0.252 ± 0.004(Gev) 



parameterized Inverse Mellin technique is used. The 
compared results will show that the CORGI approach 
indicates better consistency with available experimental 
data. An alternative "/3-fitting" method which is based 
on using the defining parameters of the Valon model, was 
also used to confirm the calculations of the parameterized 
Inverse Mellin technique. 

These calculations can be extended to the higher 
order of standard and CORGI approaches, using the 
recent analytical calculations for Wilson coefficients and 
anomalous dimensions. 



calculations. The results are tabulated in Table I. 

One way to check the validity of calculations is to ex- 
tract the quark densities using the phenomenological form 
in Eq. I|26l) . For this purpose, we went back to the /3-fitting 
method and repeated those calculations, but this time 
with the extracted values of Qq and ^jjg (a and (3 taken 
from RefQ. As a consequence, new values for parameters 
contained in Ea. i|26[l will be obtained. In order to get the 
best fit, we again multiply Ea. (|26[) by the auxiliary term 
Ax B (I — x) c ', as we did in the inverse Mellin Technique. 
We refer to this result as the "improved Inverse Mellin 
(IM) technique" . These new results in two different stan- 
dard and CORGI approaches and their comparison with 
experimental data are plotted in Fig.l and Fig. 2. As can 
be seen from Fig. 2 the result in CORGI approach will 
show better behavior in small x values. 



5. Conclusions 

Constituent quark model or valon model as a good 
candidate to describe deep inelastic scattering and to ex- 
tract sea quark densities inside the nucleon is used. The 
model bridges the gap between the bound state problem 
and the scattering problem for hadrons. In contrast to 
this model, people usually use the GLAP equations to 
evolve the parton distributions from an initial value Qq. 
The valon model which was first introduced by Hwa 
QLQI^, gives us a clear insight as to how to construct 
the hadron structure from the parton distributions. It 
has been applied to the unpolarized and polarized cases 
and gives acceptable phenomenological results. It is 
possible to use from this model in the LO and NLO 
approximation, using two different standard and CORGI 
approaches. The most important motivation for the 
CORGI approach is that by completely resumming all 
the UV logarithms one correctly generates the physical 
dependence of the moments M.(n, Q 2 ) on the DIS energy 
scale Q. 

In order to get the unknown parameters of the model 
from the fitting over the available experimental data, 
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